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GENERATION  OF  STIMULATED  BACKSCATTERED  HARMONIC  RADIATION 
FROM  INTENSE  LASER  INTERACTIONS  WITH  BEAMS  AND  PLASMAS 

I.  Introduction 

Recent  technological  advances  have  made  possible  compact  teravatt 

18  2 

laser  systems  having  high  intensities  (>  10  V/cm  ),  modest  energies 

(>  1  J),  and  short  pulses  (<  1  ps).  These  high  laser  intensities  lead 

to  a  number  of  new  laser-plasma  and  laser-electron  beam  interaction 
3-11 

phenomena.  This  paper  discusses  one  such  phenomenon,  stimulated 

backscattered  harmonic  (SBH)  radiation  generated  by  the  interaction  of 

intense,  linearly  polarized  laser  fields  with  electron  beams  or  plasmas^^ 

[see  Fig.  1].  The  intense  laser  backseat ter ing  mechanism  is  essentially 

stimulated  coherent  scattering  in  the  strong-pump  regime.  For  sufficiently 

intense  incident  laser  fields,  the  electron  quiver  velocity  becomes  highly 

relativistic.  The  high  laser  intensity,  along  with  the  induced  nonlinear 

relativistic  electron  motion,  results  in  the  generation  of  stimulated 

backscattered  radiation  at  odd  harmonics.  In  the  interaction  of  an  intense 

laser  pulse  with  a  counterstreaming  electron  beam,  SBH  radiation  is 

12-22 

generated  via  a  free  electron  laser  (FEL)  mechanism,  resulting  in  a 

13-17 

relativistic  doppler  frequency  upshift.  Hence,  a  laser-pumped  FEL 

(LPFEL)  may  utilize  both  the  harmonic  upshift  as  well  as  the  doppler 

upshift  to  generate  short  wavelength  radiation.  In  the  interaction  of  an 

intense  laser  pulse  with  a  stationary  plasma,  SBH  radiation  is  generated 

via  a  nonlinear  Raman  backseat  ter  mechanism  ’  at  odd  multiples  of  the 

fundamental  laser  frequency.  Previous  analyses  of  stimulated 

backseat ter ing  of  intense  lasers  from  electron  beams  (i.e.,  LPFELs)  or 

plasmas,  due  to  the  complex  dynamics  of  the  laser-electron  interaction, 

13-17  23 

have  been  limited  to  studies  of  the  fundamental  backscattered  mode.  ’ 

In  the  following,  a  fully  nonlinear  analysis  of  SBH  generation  is 
presented  which  is  valid  for  arbitrarily  high  pump  laser  intensities.  The 
growth  rates,  saturation  levels  (efficiencies),  and  thermal  limitations  are 
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obtained  for  the  SBH  radiation  generated  from  either  an  electron  beam  or  a 

plasma.  Differences  between  the  SBH  radiation  generated  from  an  electron 

beam  and  from  a  plasma  are  discussed.  Sufficiently  intense  pump  lasers  and 

sufficiently  cold  axial  electron  distributions  are  required  for  significant 

radiation  generation  at  high  harmonics.  The  SBH  mechanism  may  provide  a 

practical  method  for  producing  coherent  radiation  in  the  XUV  regime. 

In  conventional  FELs,  a  linearly  polarized  (planar)  static  magnetic 

viggler  may  be  used  to  generate  stimulated  harmonic  radiation  at  odd 

multiples  of  the  relativistic  doppler  upshifted  fundamental  PEL 
18-21 

frequency.  ”  In  the  LPFEL,  the  static  periodic  magnetic  viggler  is 

13-17 

replaced  by  a  counterstreaming  intense  laser  field.  For  a  given 

electron  beam  energy,  the  LPFEL  can  lead  to  substantially  shorter 

wavelength  radiation  than  a  conventional  FEL,  since  the  pump  laser 

wavelength  (typically  ~  1  ym)  is  several  orders  of  magnitude  smaller  than 

the  conventional  viggler  period  (typically  -  1  cm).  In  the  LPFEL,  the  free 

energy  driving  the  radiation  is  available  from  both  the  electromagnetic 

energy  of  the  pump  laser  and  the  kinetic  energy  of  the  electron  beam.  In 

conventional  FELs,  free  energy  is  available  from  only  the  kinetic  energy  of 

the  electron  beam  and,  hence,  the  FEL  interaction  vanishes  in  the  limit  of 

vanishing  beam  energy.  This,  however,  is  not  the  case  in  the  LPFEL 

mechanism  and  SBH  radiation  may  still  be  generated  in  the  limit  of  zero 

electron  beam  energy.  Hence,  an  intense  laser  field  interacting  «rith  a 

stationary  plasma  may  also  generate  SBH  radiation.  This  process  may  be 

viewed  as  nonlinear  Raman  backseat tering  in  the  strong-pump  limit  (or 

23 

strongly-coupled  regime). 

Significant  radiation  generation  at  high  harmonics  in  conventional 

FELs,  requires  that  the  normalized  amplitude  of  the  magnetic  viggler  field 

20 

be  somewhat  greater  than  unity.  Similarly,  in  the  LPFEL  (as  well  as  in 
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the  nonlinear  Raman  backseat taring  mechanism),  it  is  shown  that  the  growth 
rate  for  the  higher  harmonics  becomes  significant  when  the  normalized  pump 

•y 

laser  amplitude  exceeds  unity,  a^  >  1,  where  a^  »  |e|A^/m^c  is  the 
normalized  amplitude  of  the  pump  laser  vector  potential,  A^.  The 
normalized  laser  amplitude,  a^,  is  related  to  the  power,  P^,  of  a  linearly 
polarized  laser  by 

P^CGWl  .  21.5(a  r  /X  )^, 
o‘  ^  o  o  o' 

where  r^  is  the  spot  size  of  the  Gaussian  profile,  X^  is  the  laser 

wavelength  and  the  power  is  in  units  of  GW.  Physically,  a^  >  1  implies 

that  the  electron  quiver  motion  in  the  laser  field  is  highly  relativistic. 

This  may  be  seen  from  conservation  of  canonical  transverse  momentum  which, 

in  the  ID  limit  (X^  «  r^),  states  that  a^  =  where  y  is  the 

relativistic  mass  factor  and  =  Vj^/c  is  the  electron  quiver  velocity. 

The  nonlinearities  associated  with  the  highly  relativistic  electron  motion 

result  in  the  generation  of  harmonic  radiation.  In  terms  of  the  laser 

2 

intensity  (I^  *  2P^/iir^),  the  quantity  a^  is  given  by 

a  =  0.85  X  10"^X  lwmlI^^^(W/cm^], 
o  o‘  *  o  ‘ 

2 

where  X  is  in  units  of  um  and  I  is  in  units  of  V/cm  .  Relativistic 
o  o 
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electron  motion  (a^  >  1)  requires  laser  intensities  greater  than  10  W/cm 

for  wavelengths  of  ~  1  pm.  Such  intensities  are  now  available  from  compact 

laser  systems  using  a  technique  referred  to  as  chirped  pulse  amplification 

(CPA).  The  CPA  technique  allows  ultra-short  pulses  (<  1  psec)  to  be 

efficiently  amplified  in  solid  state  media  such  as  Nd:Glass,  Alexandrite 

and  TizSapphire.^  Alternatively,  KrF  excimer  laser  systems  may  be  used  to 

2 

generate  ultra-high  intensity,  subpicosecond  laser  pulses. 
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The  remainder  of  this  paper  is  organized  as  follows.  In  Section  II, 
the  generation  of  coherent  SBH  radiation  is  analyzed  using  a  relativistic, 
cold  fluid  theory.  This  theory  includes  the  nonlinear  effects  of  the  pump 
laser  to  all  orders  in  the  normalized  vector  potential,  a^.  A  dispersion 
relation  for  the  SBH  radiation  is  obtained  in  the  strong-pump  regime,  in 
which  the  effects  of  the  perturbed  electrostatic  potential  are  neglected. 
Based  on  particle  trapping  arguments,  the  saturation  level  of  the 
backscattered  harmonics  is  calculated  in  Section  III.  Limitations  of  the 
theory  due  to  thermal,  space-charge  and  collisional  effects  are  discussed 
in  Section  IV.  Examples  of  SBH  generation  using  an  electron  beam  and  a 
plasma  are  given  in  Section  V.  The  growth  rates,  saturation  efficiencies 
and  thermal  requirements  for  SBH  generation  from  electron  beams  and  plasmas 
are  summarized  in  Table  I.  This  paper  concludes  with  a  discussion  in 
Section  VI. 
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II.  Intense  Laser-Electron  Interaction  and  Stimulated  Harmonics 


The  ID  fields  associated  with  the  pump  laser,  backscattered  radiation 
and  plasma  response  are  described  by  the  transverse  vector  potential,  A, 
and  the  scalar  potential,  In  vhat  follovs  the  Coulomb  gauge  is  used, 
i.e.,  7’A  K  0,  which  in  ID  implies  A^  >  0,  where  z  is  along  the  axis  of 
propagation  of  the  SBH  radiation.  The  polarization  of  the  laser  field  is 
arbitrary.  The  normalized  potentials  satisfy  the  field  equations 


(la) 


(lb) 


222  1/2 
where  a(z,t)  =  |e|A/m^c  ,  ^(z,t)  =  |e|*/m^c  ,  Wp  =  ckp  =  (4ii|e|  n^/m^) 

is  the  ambient  plasma  frequency,  n(z,t)  is  the  plasma  density,  n^  is  the 

ambient  plasma  density,  g  «  v/c  is  the  normalized  plasma  fluid  velocity  and 

V  «  (1  -  0^  -  -  (1  +  a^)^'^^/(l  -  is  the  relativistic  factor. 

In  obtaining  the  right  side  of  Eq.  (la)  use  is  made  of  the  fact  that  the 

transverse  canonical  momentum  is  invariant  and  that  prior  to  the  laser 

pulse  interaction  the  particle  distribution  is  assumed  to  have  no 

transverse  velocity,  i.e.,  gj^  =  a/r* 

The  fluid  quantities  n,  3^,  and  y  are  assumed  to  satisfy  the  cold, 

relativistic  fluid  equations  which  can  be  written  in  the  form 


3t 


+  c 


dz 


(nP,) 


0, 


(2a) 


h  <^Pz> 


c_ 

2y  3z 


2 


a 


+  c 


31 

3z’ 


±L  1_  L. 

dt  ”  2y  3t 


2 

a 


,  ii 

z  3z’ 


(2b) 

(2c) 
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where  d/dt  ■  3/3t  +  c3  3/3z.  The  first  term  on  the  right  side  of  Eq.  (2b), 

z 

2 

proportional  to  a  ,  represents  the  ponderomotive  force.  In  the  present 
model,  thermal  effects  have  been  neglected.  This  is  valid  provided  (i)  the 
electron  quiver  velocity  is  much  greater  than  the  electron  thermal  velocity 
and  (ii)  the  thermal  energy  spread  is  sufficiently  small  so  that  electron 
trapping  in  the  plasma  wave  does  not  take  place.  The  effect  of  a  thermal 
axial  velocity  spread  on  the  wave-particle  resonance  is  discussed  in 
Section  IV.  Also,  the  ions  are  assumed  to  be  stationary. 

It  proves  convenient  to  replace  the  independent  variables  (z,t)  with 
the  independent  variables  where  h  =  z  +  ct  and  ^  =  z  -  ct.  To 

transform  Eqs.  (1)  and  (2)  from  z,t  to  ^  variables,  note  that  3/3z  =  3/3^ 
+  3/3ti  and  3/3t  =  -c(3/3^  -  3/3h).  In  the  new  variables,  Eqs.  (1)  and  (2) 
become 


iL  .  1  k2  3  ^  0 

Un  3^  A  *^p  ynj  ^ 


(3a) 


(3b) 

(3c) 


(3d) 


(3e) 
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Introducing  the  new  fluid  quantities  g  *  y(1  -  P_)»  h  ■  y(1  +  P  )  and 

z  z 

P  «  (YQn)/(Yn^),  where  y^  is  a  constant  equal  to  the  initial  relativistic 
factor  of  the  electrons  prior  to  the  laser  interaction,  Eqs.  (3)  become 


(4a) 


(4b) 


-  l^pg)  .  0. 


(4c) 


-  - 

an  " 


a  2  .fa  a  'i . 
a^  ®  ^ian  au*’ 

3  2  fa  a  'I . 

*  Hit* 


(4d) 

(4e) 


where  gh  =  1  +  a^,  y  =  (h  +  g)/2  and  =  (h  -  g)/(h  +  g). 

To  proceed  with  the  analysis,  the  radiation  is  assumed  to  consist  of  a 
large  amplitude  incident  (pump)  field  traveling  towards  the  left  and  a 
small  amplitude  backseat tered  field  traveling  towards  the  right.  The  total 
radiation  field  is  denoted  by 


a  =  a^°^n)  +  a^^^{.,n). 


where  is  the  incident  pump  (equilibrium)  field,  a^^^  is  the 

backscattered  (perturbed)  field  and  |a^°^  |  »  The  pump  field, 

a^°\  is  taken  to  be  a  function  only  of  n  (group  and  phase  velocity  equal 
to  c,  assuming  «p  «  w^),  which  implies  that  the  envelope  of  a^°^  is  non¬ 
evolving,  i.e.,  pump  depletion  effects  are  neglected.  Since  the 
backscattered  field  may  be  temporally  growing  it  is  a  function  of  both  I, 
and  q.  The  fluid  quantities  all  have  the  general  form  Q(q,  t)  =  Q^°\q)  + 

q),  where  represents  the  fluid  quantity  in  the  presence  of  only 
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the  incident  field,  is  the  perturbed  contribution  due  to  the 

backscattered  field  and 


A.  Equilibrium 

Since  the  equilibrium  quantities  are  functions  of  only  the  variable 
h,  Eqs.  (4)  become 


.  0, 

In  -  0. 


(5a) 

(5b) 

(5c) 

(5d) 


where  the  superscript  (o)  refers  to  quantities  in  the  presence  of  only  the 
incident  field,  a^°^(n),  and  =  (1  +  (a^°^)^)/h^°^ 

The  equilibrium  state  of  the  electron  fluid  in  the  presence  of  the 
pump  field  can  be  obtained  from  Eqs.  (5)  which,  in  terms  of  give 


h^°^  =  Y^d  +  Pq)  +  (6a) 

=  Y^d  +  eo)/h^°^n),  (6b) 

g^°>  =  (l  +  (a^°^n))^]/h^°^n),  (6c) 

where  the  pump  field,  a'°^(ri),  is  assumed  to  be  known,  y^  =  (1  -  3^)”^  * 
and  is  the  axial  velocity  of  the  fluid  (electron  beam  or  plasma)  in  the 
absence  of  the  pump  field.  The  self-consistent  scalar  potential  is  given 
by 
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(7a) 


where  and  C  =  kp/YQ^^^(l+0^) .  In  the  limit  that 

I  «  1»  Eq.  (7a)  reduces  to 


The  equilibrium  between  the  pump  laser  and  the  electrons  is  given  by 
the  expressions  in  Eqs.  (6)  together  with  the  solution  of  Eq.  (7a)  for  the 
potential  The  solution  of  Eq.  (7a)  is,  in  general,  highly  nonlinear. 

However,  simple  solutions  can  be  obtained  under  conditions  which  are 
relevant  to  a  wide  range  of  applications.  The  pump  laser/electron  fluid 
equilibrium  will  be  considered  in  two  limits.  In  these  limits  the 


characteristic  temporal  variation  of  the  pump  laser  envelope, 

(typically  the  pump  laser  rise  time)  is  compared  to  an  effective  plasma 
period,  (cCp)*^. 

In  the  short  pulse  limit,  x.  «  (c(c  )~^,  the  magnitude  of  is  much 

L  p 

less  than  unity,  «  1,  provided^  |a  I  <  2/(cx,lc  ),  where  a  is  the 

normalized  amplitude  of  the  pump  laser  envelope,  e.g.,  a^°^  =  a^cosk^q.  In 
this  limit,  can  be  neglected  in  Eqs.  (6). 

In  the  long  pulse  limit,  x,  »  (c£  )”^,  the  left  side  of  Eq.  (7a)  can 

L.  p 

be  neglected.  Furthermore,  it  can  be  shown  that  the  fast  oscillatory  part 
of  (on  the  laser  frequency  time  scale)  can  be  neglected^  and,  hence, 


=  (1  +  a^/2)^^^-l. 


For  applications  which  utilize  intense  pump  lasers  with  pulse  lengths 


on  the  order  of  a  picosecond,  the  short  pulse  limit  is  relevant  to 

3  16  —3 

interactions  with  electron  beams  with  densities  n  /y  «  10  cm  ,  whereas 

o  o 
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the  long  pulse  limit  is  relevant  to  interactions  with  stationary  (y^  =  1) 

plasmas  with  densities  n^  »  10  cm~  .  It  proves  convenient  to  define  the 

equilibrium  parameters  h  and  p  ,  such  that  h^°^  =  h  and  «  p  .  The 

0  0  o  o 

values  of  the  parameters  h^  and  p^  depend  on  which  of  the  above  two 
equilibriums  is  being  examined,  i.e.. 


h 


o 


-  Po>» 


e-beam  (short  pulse), 
plasma  (long  pulse). 


(8a) 


e-beam  (short  pulse), 
plasma  (long  pulse), 


(8b) 


where  =  (1  + 

Y^°^  and  n^°^  are 

2  1/2 

a^/2)  .  Hence,  the  equilibrium  fluid 

given  by 

quantities 

«<»>  .  [ 

^  1  *  r'n. 

(9a) 

-  [ 

(9b) 

n*”)  .  ( 

"<.v2Vo)[''o  *  1 

(9c) 

Physically,  the  difference  between  the  plasma  (high  density)  equilibrium 
and  the  electron  beam  (low  density)  equilibrium  is  due  to  the  fact  that  in 
the  high  density  regime  the  ponderomotive  force  associated  with  the  laser 
envelope  variation  is  balanced  by  the  space-charge  force,  whereas  in  the 
low  density  regime  the  ponderomotive  force  dominates. 

B.  Backscattered  Radiation 

To  analyze  the  SBH  radiation,  Eqs.  (4)  are  expanded  about  the 

equilibrium  state  given  in  Eqs.  (6).  In  the  following,  the  perturbed 

12  20  23 

electrostatic  field,  is  neglected.  This  is  valid  provided  ’  ’ 
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the  temporal  growth  rate  of  the  backscattered  harmonic  radiation  is  much 

3/2 

greater  than  the  relativistically  corrected  plasma  frequency,  •  In 

this  strong-pump  regime,  the  SBH  radiation  is  completely  described  by  the 
fluid  equations 


P  k  -  Id''  -  -  h  It 

P  In  -  It]'  -  -  '  In 


(10a) 


(10b) 


along  with  the  wave  equation,  Eq.  (4a),  and  the  equilibrium  given  by  Eq. 

(8). 

The  self-consistent  closed  set  of  perturbed  quantities,  h^^\  and 

a^^\  satisfy 


(In  -  """Id-' 


k^p  ^ 
p  o ' 


(1) 

_  _  2_  i.  a<^>) 

(11a) 

o 

_  ^  L.  h(l), 

(11b) 

(1) 

_ 

(11c) 

3n3^  4y, 


where  w(t|)  =  (1  +  (a^°^)^)/h^.  To  proceed,  the  perturbed  quantities  are 
represented  by  the  form  =  Q2^(h)e^^^/2  +  c.c.,  where  k  is  complex  and 

the  amplitude  is  a  function  of  q.  Using  this  representation  the  perturbed 
amplitudes  are  given  by 


(1^  -  iku(n))hi(n)  =  -  ^  a^°^ti)*aj(n), 

(In  -  "-"C'd-'i"’'  ■  -  r  "o  In"!""’ 


(12a) 


(12b) 
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(12c) 


Ih  ■  '  4?" 

o  0 


Equations  (12a)  and  (12b)  have  the  solutions 


h^Ol)  =  e  ^  a^°)(n' )-aj(n' ))dn' ,  (13a) 


Pi  (11) 


^  ^  h^(n)  *  e*®*"*/  .-‘®<''')  (-  w(n')  hj(,-)]dn', 


(13b) 


•  I 

where  0(ri)  =  kj  u(ri")dh”. 


Taking  the  incident  field  to  be  a' 


a  cosk  h  and  the  complex 
~o  o 


amplitude  of  the  backseat tered  harmonic  field  to  be  a^  >  a^exp(iAkh),  where 
Ak  is  complex,  the  solutions  to  Eqs.  (13a)  and  (13b)  become 

.  V  -  -'e.l*'”] 

1  h„  *o  *iZj  fk  -  k  (l*2e)  -  4k] 

t,n=-®  '  ° 


>exp{i[k^(l  +  2(e+n))  +  AkJn), 


(lAa) 


Pl(>l) 


2kk  A  V  .  ... 

=  — T  io'  ^1  /  7= - 

h^  ^  ^  ^  "  k  -  k  (U 

o  4  _  _  (o' 


e,n=-® 


h(b)  -  Jt.l<b)) 

(ic  -  k^(l+2e)  -  Ak]^ 


'exp{i[k^(l+2(e+n))  +  Akjn, 


(14b) 


—  2  2  2  2 

where  k  =  obtaining  the  expressions  for 


h.  ()i)  and  Pi()l),  the  Bessel  identity,  exp(ibsinx)  »  u  (b)exp(inx),  was 
used.  Also,  near  resonance  the  first  term  on  the  right  of  Eq.  (13b)  has 
been  neglected  and  the  initial  value  terms  in  Eqs.  (13a)  and  (13b)  have 
been  neglected  compared  to  the  exponentially  growing  modes. 


C.  Dispersion  Relation 

The  dispersion  relation  for  the  SBH  modes  may  be  obtained  from  the 
perturbed  wave  equation.  Substituting  Eq.  (14b)  into  Eq.  (12c)  and 
requiring  that  both  sides  of  the  equation  have  the  same  q  dependence, 
yields 


9 

4kAk  + 

P 


p  k^kk 
o  P 

rh2 
o  o 


(•'e  -  ■'e.i)' 


[k  -  (2e+l)k^-Ak]* 


(15) 


where  ^  =  ®o^l  been  assumed.  The  frequency,  u,  and  temporal 

growth  rate,  T,  may  be  determined  from  the  above  expression  by  setting 
k  =  (o/c  -  Ak.  The  backseat tered  radiation  has  the  phase  dependence  a^^^  ~ 
exp[i(w/c)^  +  i2ctAk],  where  the  temporal  growth  rate  is  given  by  T  =  - 
2cIm(Z^).  The  dispersion  relation  indicates  that  the  resonant  frequency  of 
the  mode  is 


(I)  -  N  M  <«>  , 
o  o’ 


(16) 


where  N  =  2e.+l  is  the  harmonic  number  and  M  =  (h  /y,  )^  is  the  frequency 

o  O  10 

multiplication  factor  which  is  dependent  on  the  particular  equilibrium 
being  examined,  i.e.. 


M  = 
o 


I'Y^d  +  +  aJ/Z),  e-beam, 

il,  plasma. 


(17) 
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Assuming  |«/c|  »  |Ak|  »  »  the  dispersion  relation  reduces  to 


.3  »oWo^<\> 

OK  *  —  9  > 

Y  (1  +  M 
'o'  o' 


(18) 


2  20 
where  =  b^[J^(b^  ^  ■  *^£+1  (b^)]  is  the  harmonic  coupling  function  and 

2  2 

b^  -  (2£  +  The  temporal  growth  rate  is 


r  =  /Sc 


P  k  k  M 
o  p  o  o  £ 


1/3 


(19a) 


A  plot  of  the  function  F„  versus  b  =  a  /4yT  for  N  =  l+2£  =  1,3,5 _ 19 

£  o  o  ±o 

is  shown  in  Fig.  2.  For  a  relativistic  electron  beam  in  which  h^  =  » 


Yj^q,  the  growth  rate  is  given  by 


1/3 


r./3c[kWp,/4r5]  . 


(19b) 


For  a  plasma,  the  growth  rate  is  given  by 


r  =  /3c[k2k„Fj/4r,„J 


.1/3 


(19c) 


The  asymptotic  value  of  the  harmonic  coupling  function  F^(b^)  for 

2 

a^  »  1  and  £  >>  1  can  easily  be  found.  In  this  asymptotic  limit,  b^  =  £ 

f  -2/3 

and  J^(£)  -  J^^i(^)  =  J^(t)  =  0.41U  Hence, 

F^  =  0.169/£^''^.  (20a) 

2  -1/9 

Notice  that  asymptotically  (£  »  1,  a^  »  1),  T  -  £  ,  i.e.,  the  growth 

rate  is  a  weakly  decreasing  function  of  £. 

For  small  arguments,  b^  <  1,  the  harmonic  coupling  function  is  given 
by 


14 


(20b) 


,  2 

Hence,  in  the  limit  (2t+l)a^/2  «  1,  the  growth  rate  is  a  rapidly 

2 

decreasing  function  of  £.  Notice  that  for  aV2  «  1  and  t  -  0,  Eqs.  (19c) 

and  (20c)  indicate  the  growth  rate  of  the  fundamental  backseat tered 

radiation  from  a  plasma  is  F  =  /3(«p«^a^/16)  ,  which  is  the  standard 

result  for  the  Raman  backscatter  instability  in  the  strong-pump  or 

2 

strongly-coupled  regime  for  a  linearly  polarized,  low-power  (a^  «  1) 

,  23 

laser. 
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III.  Saturation 

The  present  analysis  indicates  that  the  SBH  radiation  grows 
exponentially,  with  a  growth  rate  given  by  Eq.  (19).  This  growth  continues 
until  nonlinear  effects  (i.e.,  particle  trapping)  limit  the  amplitude.  The 
saturated  value  of  the  backscattered  radiation  may  be  determined  from 
arguments  based  on  electron  trapping  in  the  ponderomotive  wave.  The  growth 
of  the  radiation  will  cease  when  the  generalized  ponderomotive  wave  traps 
the  electrons.  For  small  amplitude  pump  fields,  |a^|  «  1,  the 
ponderomotive  wave  is  proportional  to  the  product  of  the  stimulated 
backscattered  radiation  and  the  pump  field,  •  a^°^.  In  the  present 

analysis  a  generalized  ponderomotive  wave  is  obtained  for  arbitrarily  large 
values  of  |a^|.  The  nonlinear  effects  of  the  pump  field,  in  the 
generalized  ponderomotive  wave,  are  included  to  all  orders. 

The  combined  action  of  the  pump  laser  and  the  stimulated  backscattered 
radiation  results  in  a  generalized  ponderomotive  wave  which  produces  a 
fluid  wave,  with  the  same  phase, 

^  exp  |i[(2t  +  1  +  2n)k^  +  ^)]n  +  iktl,  (21) 

as  indicated  by  Eq.  (14b),  where  p^^^  =  Pj^(ri)exp(ik^)/2  +  c.c.  and  ^  are 
constant  coefficients.  The  generalized  ponderomotive  wave  consists  of  a 
sum  of  individual  modes  characterized  by  the  mode  numbers  (?.,n).  Since 
1.  =  z-ct,  q  =  z+ct  and  k  =  (o/c  -  dk,  the  normalized  phase  velocity  for  a 
particular  generalized  ponderomotive  mode  (t,n)  in  Eq.  (21)  is  given  by 

w/c  -  (21  +  2n  +  l)k  -  2dk 
=  w/c  +  (2t  +  2n  +  l)ir  • 
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At  a  particular  resonant  frequency  a  *  (21  +  of  the  mode,  the 

phase  velocity  in  Eq.  (22a)  becomes 


ep(t,n) 


(M  -  l)(2t  +  l)k  -  2nk  -  2Ak 
o  _  o _ o _ 

(M^  +  l)(2e  +  l)k^  +  2nk^ 


(22b) 


Particle  trapping  for  a  particular  n  mode  occurs  vhen  the  longitudinal 
fluid  velocity  of  that  mode  is  equal  to  the  phase  velocity  in  Eq.  (22b). 
Trapping  occurs  first  for  the  n  «  0  mode,  since  it  is  the  mode  with  phase 
velocity  closest  to  the  equilibrium  velocity  of  the  electrons. 

The  longitudinal  fluid  velocity  of  the  electrons  is  given  by  3^  = 
(h-g)/(h+g),  which  can  be  written  as 


6  -  (23) 

*  h^  +  (1  +  a^) 

In  terms  of  equilibrium  and  perturbed  quantities,  3^  =  3^°^  where 

the  average  value  of  3^°^  is  <3^°^>  -  (M  -1)/(M  +1).  Expanding  Eq.  (23), 
the  leading  order  contribution  to  the  perturbed  velocity  at  resonance  is 
given  by 


»(1)  - 


4H 


,(1) 


‘  "o) 


2  h 


(24) 


where  h^^^  is  composed  of  a  collection  of  n  modes,  the  phase  velocity  of 


the  n^^  mode  given  in  Eq.  (22b).  Using  Eq.  (14a),  the  perturbed  fluid 


velocity  of  the  n^^  mode  is  given  by 


>0<1  » 


(bF^) 


(25) 


Particle  trapping  occurs  first  for  the  n  =  0  mode,  which  is  the  mode 
with  phase  velocity  closest  to  the  equilibrium  fluid  velocity, 
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(26a) 


e„(e,n  -  0)  .  +  SB, 

P  *  P 

where 


63p  =  -  2Ak/(2£  +  1)(1  + 


(26b) 


The  saturated  level  of  the  radiation  field  may  be  obtained  from  the 

condition  for  deeply  trapped  electrons,  |  =  2|&3  |,  which  gives  for 

z  p 

the  n  <=  0  mode, 


4^j(2e  +  ml 

a  (1  +  M  )^ 
o'  o' 


(bF^) 


(27) 


A  numerical  factor  of  order  unity  has  been  neglected  on  the  left  side  of 
Eq.  (27).  Using  the  dispersion  relation,  the  ratio  of  the  radiation  power 
in  the  harmonic  at  saturation  to  that  in  the  pump  laser,  Pjj/Fq  * 
la^l^/a^,  is  given  by 


p  k‘(l  +  M  ) 
*^o  p' _  o_ 


,  2„l/2 
Y  k  M 
o  o  o 


4/3  U/3 


(28) 


1/3  2  2 

A  plot  of  the  function  /b^  versus  b^  =  for  N  =  l+2e  = 

1,3, 5... 19  is  shown  in  Fig.  3. 

2  -10/9  2 

Note  that  asymptotically  (t  »  1,  a^  »  1),  Fj^/P^  -  I  /J^. 

Furthermore,  the  efficiency  may  be  substantially  increased  by  operating  in 

2  2, 

a  regime  in  which  b^  =  (1  +  2t)a^/4Yj^Q  is  near  a  zero  of  J^,  i.e.,  J^(b^)  = 
0.  Physically,  the  regime  for  which  J^(b^)  =  0  for  a  particular  C  mode 
corresponds  to  minimizing  the  perturbed  fluid  velocity  n  =  0)  of 

that  mode.  This  implies  that  a  larger  amplitude  of  backscattered  radiation 
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will  be  reached  before  saturation  occurs.  At  J^(bj^)  -  0,  Eq.  (28)  no 

longer  applies,  and  saturation  due  to  particle  trapping  must  be  determined 

by  consideration  of  the  higher  n  modes.  In  general,  Eq.  (28)  is  valid 

provided  <1.  If  this  is  violated,  then  saturation  will  be  caused  by 

some  other  mechanism  besides  particle  trapping,  i.e.,  plasma  thermalization 

or  pump  depletion.  In  particular,  notice  that  for  a  plasma  with  e  «  0  and 

a^  «  1,  Eq.  (28)  indicates  Pjj/Pq  ®  hence,  Pjj/Pq  <  1 

implies  a  >  k^/2k^. 

o  p  o 

It  is  also  convenient  to  define  an  electronic  efficiency,  which  is 
the  ratio  of  the  backscattered  radiation  power  to  the  electron  beam  power, 


u2„2,  2  lA  i2 

H  N  k  a, 
o  0*1' 

2k^e  (y  -  1) 
p  o'  o  ' 


(29a) 


Using  Eq.  (28),  the  electronic  efficiency  may  be  written  as 


8J^M  p  (1  -  1/y  )  Nk  ’ 
O  O  '  '  0  ^  o 


(29b) 


where  |Ak|  =  r//3c  is  given  by  Eq.  (19).  Notice  that  for  the  fundamental 

2  _ 

(£,  =  0)  in  the  limit  y^  »  1  and  a^  «  1,  hg  =  r/2/3(«)^,  as  is  the  case  for 

12  —  2  2  3  1/3 

conventional  strong-pump  FELs,  where  F  =  /3(«p(i)^ap/16Yg) 
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IV.  Thermal,  Space-Charge  and  Collisional  Effects 


A.  Thermal  Effects 

The  above  results  were  obtained  using  cold  fluid  theory,  i.e., 
electron  thermal  effects  were  neglected.  For  a  sufficiently  thermal 
electron  distribution,  however,  it  is  possible  for  the  longitudinal  energy 
spread  to  become  large  enough  so  as  to  degrade  the  resonant  interaction 
between  the  backscattered  wave  and  the  electron  distribution.  This  thermal 
electron  interaction  regime  corresponds  to  a  weak  resonant  instability  in 
which  the  growth  rate  of  the  backscattered  radiation  is  greatly  reduced. 

It  is  possible  to  estimate  how  large  a  thermal  electron  velocity 
spread  can  be  tolerated  before  the  resonant  wave-particle  interaction  is 
degraded.  The  growth  rate  of  the  cold  electron  instability  is  determined 
largely  by  the  resonant  denominator,  D,  in  the  dispersion  relation,  Eq. 
(15), 

D  *  (k  -  (2e+l)kQ  -  Akl^,  (30) 

where  k  =  k/M^  and  is  given  in  Eq.  (17).  The  effects  of  a  longitudinal 
thermal  velocity  spread  on  the  resonance  may  be  estimated  by  letting 
-»  +  p^j^  in  the  expression  for  D,  where  ^  resonance, 

(o/c  =  M^(2C.+l)k^,  one  finds 

D  =  I2ro(l+2e)k^P^jj  +  (Ul/M^)Akl2.  (31) 


Hence,  in  order  to  neglect  for  the  effects  of  the  thermal  velocity  spread 
P^j^,  it  is  necessary  that 


(UM^)|Ak|/k^ 

2yh\^{2Ul) 


(32) 
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Using  the  dispersion  relation,  Eq.  (32)  yields  the  condition 


^th  2 

2y  N 
o 


o  o  o 


1/3 


(33) 


1  /T  9  9 

A  plot  of  the  function  F]:  /N  versus  b  *  a  /4y7  for  N  =  l+2£  =  1,3, 5.  ..19 

£  0  0  10 

2 

is  shown  in  Fig.  4.  Asymptotically  (t  »  1  and  a^  »  1),  the  right  side  of 

Eq.  (33)  scales  as  For  a  stationary  plasma,  the  thermal  energy  of 

2  2 

the  plasma  is  AE^^  =  m^c  0^^/2.  For  an  electron  beam  with  initial  velocity 


'th 


0Q  »  and  Yq  »  1»  the  normalized  energy  spread  is  given  by  Ay^Jj/y^  = 

2  2 
Yo^th*  Note  that  for  the  fundamental  (£  =  0)  in  the  limit  y^^  »  1  and  a^ 

«  1,  the  thermal  requirement  is  Ay^j^/y^  «  as  is  the  case  for 

12 

conventional  strong-pump  FELs  operating  at  the  fundamental.  The  usual 
requirement  regarding  FELs  operating  at  the  fundamental,  Ay^j^/y^  « 
however,  does  not  apply  in  general  to  harmonic  generation. 


B.  Space-Charge  Effects 

The  expressions  for  the  growth  rate  of  the  SBH  radiation,  the 
amplitude  of  the  radiation  at  saturation  and  the  allowable  thermal  energy 
spread  are  all  increasing  functions  of  the  electron  density.  Hence, 
optimal  generation  of  the  SBH  radiation  implies  operating  in  a  regime  of 
high  electron  density.  However,  the  fluid  theory  of  Section  II  used  to 
describe  the  backscattered  radiation  assumed  the  strong-pump  limit,  i.e., 
assumed  that  the  electron  density  was  sufficiently  low  so  that  the  effects 
of  the  perturbed  electrostatic  potential  may  be  neglected.  This  is  valid 
provided  the  relativistic  plasma  frequency  is  small  compared  to  the  growth 
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rate  of  the  radiation.  More  specifically,  it  can  be  shown  that  the  strong- 
pump  limit  is  valid  provided 


P  k  M  . 

°  ,  «  |dk|^.  (34) 

Y  (1+M  ) 

'o'  o' 

This  inequality  gives  an  upper  limit  to  the  electron  density  for  which  the 
strong-pump  limit  remains  valid.  Using  the  dispersion  relation,  Eq.  (18), 
gives 


k^/k 

P 


2 

o 


«  Y  (1+M  )^F?/p  M  . 
o  o  £  o  o 


(35) 


2 

Note  that  asymptotically  (£  »  1  and  a^  »  1),  the  right  side  of  Eq.  (35) 
-2/3 

scales  as  i  .  Strictly  speaking,  the  results  of  the  strong-pump  theory 
obtained  in  the  previous  sections  are  valid  provided  Eq.  (35)  is  satisfied. 


C.  Collisional  Effects 

The  effects  of  collisions  on  the  backscattered  harmonics  may  be 

neglected  provided  the  damping  rate  of  the  radiation  field  due  to 
3  2  2 

collisions,  v  =  («p/«  )'Vgj,  is  small  compared  to  the  growth  rate  of  the 

radiation  field,  F,  where  is  the  electron-ion  collision  frequency, 

2  2  24 

i.e.,  r  »  (Wp/d)  The  relativistic  collision  frequency  is  given  by 

-3  -2 

v^jlsec”^]  =  4.3  X  10~^^n^Icm'^J(lnA)  j  »  ^^6) 


where  InA  is  the  Coulomb  logarithm.  The  collision  frequency  will  be 
largest  for  electrons  in  a  dense  plasma.  For  plasma  electrons  undergoing 
relativistic  quiver  motion  in  the  intense  pump  laser  field, 


^lo 

1  ^  1 

f  2/,  2  ) 
a  / 4y, 

1  0  'ioj 

IcosZk  h 
' _ o' 

[cos^c 

+  1 

1  2.,. 
cos  2k 
'  c 

n] 

3/2’ 


(37) 
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An  upper  limit  on  Eq.  (37)  is  given  by 


(38) 


2  2 
In  the  limit  3^/2  »  1,  the  right  side  of  Eq.  (38)  reduces  to  24/a^. 

2 

Hence,  for  intense  laser  fields  where  a^  >  1,  the  condition  that  F  v  is 
easily  satisfied. 
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V.  Examples 


The  SBH  instability  described  in  the  previous  sections  may  be  used  as 

a  mechanism  to  generate  coherent  radiation  in  the  XUV  regime. 

Concept ionally,  a  device  may  be  designed  to  amplify  a  small  XUV  input 

signal  injected  into  a  plasma  or  a  co-propagating  electron  beam.  The  input 

XUV  signal  may  be  obtained  by  the  output  of  an  incoherent  source,  such  as  a 

flash  lamp,  or  by  using  the  incoherent  single  particle  (spontaneous) 

radiation  generated  by  the  interaction  of  the  pump  laser  with  the  electron 

distribution.  The  XUV  radiation  will  be  amplified  as  it  propagates,  via 

the  interaction  with  the  intense,  counter-propagating  pump  laser  field. 

Two  examples  will  be  discussed;  one  utilizing  a  stationary  plasma  and  the 

other  utilizing  a  relativistic  electron  beam.  In  both  examples,  the  pump 

laser  is  taken  to  be  an  intense  Nd:Glass  laser  (X  =1  pm)  with  an 

18  2 

intensity  of  =  9.7  x  10  W/cm  (a  =  2.6)  and  a  spot  size  of  r  = 
o  o  o 

10  pm,  which  gives  a  pump  laser  power  of  =  15  TV. 

A.  Backseat tering  from  a  Plasma 

Consider  the  amplification  of  coherent  radiation  using  a  plasma  of 
19  -3 

density  n^  =10  cm  ,  which  corresponds  to  fully  ionized  H2  at  about 

5-8 

0.1  atm.  Neglecting  the  effects  of  relativistic  optical  guiding,  ~  the 

interaction  length  is  approximately  either  two  vacuum  Rayleigh  lengths,  2Zj^ 
2 

=  2nr^/X^  =  630  pm,  or  one-half  the  pump  laser  pulse  length  (which  is  150 
pm  for  a  1  psec  laser  pulse),  whichever  is  shorter.  As  an  example, 
consider  the  amplification  of  the  third,  f  =  1,  harmonic  (and  the  fifth, 

O  0 

t  =  2,  harmonic)  at  a  wavelength  of  X  =  3300  A  (2000  A).  The  e-folding 

length  is  c/T  =  1.8  pm  (2.0  pm).  At  saturation,  the  ratio  of  the  harmonic 

-A  -4 

power  to  the  pump  laser  power  is  =  1*0  ^  10  (3.8  x  10  ),  which 

implies  a  saturation  power  of  Pj^  =  1.5  GU  (5.8  GW).  (The  larger  saturated 
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2 

power  for  the  fifth  harmonic  is  due  to  the  factor  1/J^.)  The  pulse  length 

of  the  SBH  radiation  is  approximately  the  transit  length  of  the  pump  laser 

5-8 

pulse  through  the  plasma.  Optical  guiding  may  substantially  increase 

this  length  beyond  the  vacuum  diffraction  limit  of  >  630  urn.  The 

thermal  requirement  on  the  longitudinal  energy  spread  is  <  77  eV 

(22  eV).  Plasmas  with  sufficiently  cold  longitudinal  temperatures  may  be 

25 

produced  by  laser-induced  ionization. 

B.  Backseat  taring  from  an  Electron  Beam 

Consider  amplification  of  coherent  radiation  using  an  intense  electron 

beam  with  a  current  of  15  A,  a  beam  radius  of  10  ym  (a  current  density  of 
2 

4.8  MA/cm  )  and  an  energy  of  250  keV  (y^  =  1.5).  The  interaction  length  is 
approximately  either  two  vacuum  Rayleigh  lengths,  22^  =  630  ym,  or  one-half 
the  laser  pulse  length,  whichever  is  shorter.  As  an  example,  consider  the 
amplification  of  the  third,  £.  =  1,  harmonic  at  a  wavelength  of  X  *  X^/3M^  = 

O 

2200  A.  The  e-folding  length  is  c/T  =  31  ym.  The  thermal  requirement  on 
the  longitudinal  energy  spread  is  br/(y^  -  1)  <  0.18Z.  The  pulse  length  of 
the  SBH  radiation  is  approximately  the  transit  length  of  the  pump  laser 
pulse  through  the  electron  beam,  i.e.,  2Zj^  =  630  ym.  At  saturation,  the 
ratio  of  the  third  harmonic  power  to  the  pump  laser  power  is  Pjj/Pq  = 

_Q 

1.1  X  10  ,  which  implies  a  saturation  power  of  Pj^  =  17  kV.  The  electronic 

efficiency  at  saturation  is  =  0.87Z. 
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VI.  Discussion 


The  generation  of  coherent,  stimulated  backscattered  harmonic  (SBH) 

radiation  by  the  interaction  of  an  intense,  linearly  polarized  laser  pulse 

with  an  electron  beam  or  a  plasma  has  been  analyzed  using  relativistic, 

cold  fluid  equations.  This  theory  includes  the  nonlinear  effects  of  the 

pump  laser  to  all  orders  in  the  normalized  vector  potential,  a^.  A 

dispersion  relation  for  the  SBH  radiation  was  obtained  in  the  strong-pump 

regime.  In  this  regime,  the  effects  of  the  perturbed  electrostatic 

potential  are  neglected.  The  resonant  frequency  of  the  mode, 

corresponding  to  the  harmonic  number  N  =  2t.+l,  is  given  by  w/c  =  NM^k^, 

where  is  the  frequency  multiplication  factor  arising  from  the 

relativistic  doppler  upshift.  The  temporal  growth  rate  is  given  by  Eq. 

(19).  Significant  high  harmonic  generation,  N  »  1,  requires  intense  pump 

lasers,  a  >  1.  The  SBH  radiation  from  an  electron  beam  may  be  viewed  as  a 
0 

12-22 

nonlinear  laser-pumped  PEL  mechanism,  whereas  the  SBH  radiation  from  a 
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plasma  may  be  viewed  as  a  nonlinear  Raman  instability  mechanism,  both  in 
the  strong-pump  regime. 

The  saturation  amplitude  of  the  SBH  radiation  was  calculated  based  on 

particle  trapping  arguments.  Particle  trapping  (wavebreaking)  occurs  when 

the  amplitude  of  the  perturbed  longitudinal  fluid  velocity,  becomes 

equal  to  the  phase  velocity,  B  .  The  perturbed  fluid  velocity, 

consists  of  a  sum  of  individual  waves  characterized  by  the  mode  numbers 

(£,,n).  For  a  particular  t  resonance,  particle  trapping  occurs  first  for 

the  n  =  0  mode,  which  is  the  component  of  with  phase  velocity  closest 

to  the  equilibrium  longitudinal  electron  velocity,  <B^^^>.  The  resulting 

saturation  amplitude  of  the  backscattered  radiation  is  given  by  Eq.  (26). 

Furthermore,  the  saturation  amplitude  of  the  backscattered  radiation  may  be 

2  2 

enhanced  by  operating  in  a  regime  in  which  b^  =  (2£.+l)a^/4Yj^^  is  close  to  a 
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zero  of  the  Bessel  function  J^(b^).  Physically,  =  0  corresponds  to 
minimizing  the  amplitude  of  the  n  «  0  component  of  the  perturbed  fluid 
velocity  When  this  occurs,  saturation  may  be  determined  from 

consideration  of  the  next  higher  n  mode.  It  should  also  be  pointed  out 
that  particle  trapping  only  leads  to  saturation  for  sufficiently  intense 
pump  laser  amplitudes  for  which  Pjj/Pq  ^  1*  When  this  inequality  is  not 
satisfied,  saturation  occurs  by  other  nonlinear  effects,  such  as  electron 
thermalization  or  pump  depletion. 

The  most  stringent  constraint  on  the  production  of  SBH  radiation  is 
the  restriction  on  the  longitudinal  thermal  velocity  spread.  If  the 
longitudinal  electron  temperature  is  sufficiently  high,  the  wave-particle 
resonance  is  degraded  and  the  growth  rate  of  the  backscattered  radiation  is 
greatly  reduced.  The  expression  for  the  allowable  thermal  velocity  spread 
is  given  by  Eq.  (33).  In  general,  this  implies  that  very  small  thermal 
spreads  are  required  in  order  to  generate  high  order  harmonic  radiation. 

The  expressions  for  the  growth  rate,  saturation  amplitude  and  allowable 
thermal  velocity  spread  indicate  that  high  electron  densities  are  required. 
This  appears  to  favor  the  use  of  stationary  plasmas  over  that  of 
relativistic  electron  beams.  However,  the  use  of  an  electron  beam  has  an 
advantage  in  that  the  frequency  of  the  harmonic  radiation  is 
relativistically  upshifted  (in  addition  to  the  harmonic  upshift),  which 
implies  the  radiation  frequency  may  be  tuned  by  adjusting  the  energy  of  the 
electron  beam  or  the  amplitude  of  the  pump  laser,  as  indicated  by  Eq.  (16). 
The  growth  rates,  saturation  efficiencies  and  thermal  spread  requirements 
for  the  SBH  radiation  are  summarized  in  Table  I.  As  the  capability  for 
producing  dense  plasmas  and  electron  beams  with  small  thermal  spreads 
improves,  along  with  future  advances  in  ultra-high  power  laser  technology. 
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SBH  generation  may  provide  a  practical  method  for  producing  coherent 
radiation  in  the  XUV  regime. 
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Table  I. 


Growth  Rates,  Efficiencies  and  Thermal  Requirements  for 
Stimulated  Backscattered  Harmonic  Generation^ 


Growth  Rate 
r/wo 


Laser  Eff.^ 

Pi/Po 


Electron.  Eff. 


Thermal  Spread 
AEth/moc^  < 


Energy  Spread 
A7tfc/7o  < 


^  u>o  is  the  pump  laser  frequency,  tjp  is  the  plasma  frequency,  70  is  the  initial  relativistic 
factor,  oo  =  |c|i4o/moc*  is  the  normalized  pump  laser  amplitude,  7j^o  =  (1  +  “o  N  = 

(2^4-1)  is  the  harmonic  number.  Ft  =  6[J/(6)  —  Ji+i{b)]^  is  the  harmonic  coupling  function, 
b  =  TV 00/47^0  and  Mq  is  the  frequency  amplification  factor,  u;  =  NMou;o. 

^  Formulas  valid  provided  oo  >  u;^/2u;q. 


Intense  Pump  Laser,  cJq 


Electron  Beam 
or 

Plasma 


Backscattered  Harmonics,  Wq 


NMqWo 


Fig.  1  Schematic  of  the  incident  pump  laser  field,  ~  exp(ia)^ri/c) ,  of 
frequency  «^,  and  the  backscattered  harmonic  radiation  field, 

~  exp(ico^c),  of  frequency  «,  as  obtained  via  stimulated 
scattering  from  an  electron  beam  or  a  plasma,  where  q  =  z  +  ct  and 
^  =  z  -  ct. 
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{allA)/{l  +  al/2) 

1/3 

The  function  (proportional  to  the  growth  rate  F)  versus  the 

2  2 

parameter  =  (a^/4)/(l  +  a^/2)  for  the  harmonics  N  =  1,3,5...!' 


The  function 

P„/P  )  versus 
N  o 


(a5/4)/(l  +  aim 
1/3 

^  /b^  (proportional  to  the  power  at  saturation 

2  2 

the  parameter  =  (ap/A)/(l  +  3^/2)  for  the 


harmonics  N  =  1,3,5. ..19. 


{allA)/{l^al/2) 

1/3 

The  function  /N  (proportional  to  the  upper  limit  of  the 

O 

thermal  velocity  spread,  versus  the  parameter  =  (aV4)/(l 

2 

+  a  /2)  for  the  harmonics  N  =  1,3,5 _ 19. 


